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We study the consequences of the large dimuon asymmetry observed at DO. 
(— i ' Physics beyond the standard model (SM) in B s — B s mixing is required to explain 

the data. We first present a detailed analysis for model independent constraints 

CO 

\ on physics beyond the SM, and then study the implications for theoretical models 

OO 



O 
O 



X 



which modify the SM results in different ways, such as Z' with FCNC and R-parity 
violating SUSY contributions. 

PACS numbers: 



I. INTRODUCTION 

The DO Collaboration, with 6.1 fb _1 of data, has recently reported evidence for an anoma- 
lously large CP violation in the like-sign dimuon charge asymmetry [l| which is attributed 
to semileptonic decays of b hadrons defined by 

A b — __2 o (I \ 

sl " N++ + N b - ' { > 

where Njj~ + (N b ) is the number of events with two b hadrons (b b) decaying semileptonically 
into /I+/I+X (jTfj,-X). The DO result fl, A b sl = -(9.57 ± 2.51 ± 1.46) x 10~ 3 with the first 
error being statistical and the second systematic, is 3.2er away from the standard model 
(SM) prediction of -0.2 x 10~ 3 Q. 
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Agj is related to the asymmetries a s { s in B d and B s decays by 

where Z 9 = 1/(1 - j/J) - 1/(1 + x 2 q ) with ^ = Aiy2r, and x q = AM q /T q , and f d and /, 
stand for the fragmentation fractions for b — > B d and b — >• _B S , respectively. Using known 



values of 



3j f d = 0.323 ± 0.037 and f s = 0.118 ± 0.015, x d = 0.774 ± 0.008 and y d = 0, 



and x s = 26.2 ± 0.5 and y s = 0.046 ± 0.027, one has 

A b sl = (0.506 ± 0.043)4 + (0.494 ± 0.043)4 . (3) 



In the above is the "wrong-charge" asymmetry, 



(4) 



81 r(5 g ^+x) + r(5 ? ^-x) 

Using the current experimental value [l], Ql afj = —0.0047 ± 0.0046 which is consistent 
with zero and also with the SM prediction [2j = (-4.8lJ;°) x 10" 4 within errors, one 
requires [l| 

a s sl = -0.0146 ± 0.0075 , (5) 

to obtain the DO value of A^. This value is much larger than the SM prediction {2]] (2.1 ± 
0.6) x 10~ 5 for a*. 

The CDF j^J measurement of A^, using only 1.6 fb" 1 of data, has a central value which 
is positive, A^ = (8.0 ± 9.0 ± 6.8) x 10 -3 , but is still compatible with the DO measurement 
at the 1.5cr level because its uncertainties are 4 times larger than those of DO. Combining 
in quadrature (including the systematic errors) the DO and CDF results for A^, one finds 
A b sl ~ -(8.5 ± 2.8) x 10~ 3 which is still 3o" away from the SM value. 

DO [6J also performed a direct measurement of a^, but the result does not show any 
deviation from the SM, although the error bars are quite large: a s sl = -(1.7±9.1^; 4 ) x 10~ 3 . 
This value is much smaller than the one given in Eq. ([5]). We could combine all these results 
to obtain an average value {a^)^ ~ — (12. 7 ± 5.0) x 10 -3 . 

Even though the inclusion of the CDF dimuon asymmetry and the DO semileptonic wrong- 
charge asymmetry reduces the deviation in derived from the DO dimuon asymmetry, the 
above result is still about 2.5cr away from the SM value [2]] of (g^Jsm- If confirmed, it is an 



indication of new physics 



data have been made 



jeyond the SM [?-l2l|. Several theoretical attempts to explain the 



3J-|2l| 



We note that there could be charm contamination. It is known that there is D° — D° 
mixing, and this will modify the asymmetry measured by DO unless D mesons are completely 
eliminated as a possible source of muons in their data sample. To take this contamination 
into account one would add terms proportional to f°Z°, related to charm contribution, to 
the formula in Eq. 02]). Here is analogous to Zj and is determined by the D° — D° mixing 
parameters x u and jy„, and is the fraction of direct D° and D° production in pp collisions. 



Using current values 



22| of x u = (0.98±g;it)% and y u = (0.83 ± 0.16)% for D° - D° mixing, 



the factor can be determined. With the central values we obtain Z^/Z^ m 7.2, which is 
a large number. In order to obtain a reliable result, the parameter f° should be carefully 
evaluated for the muon selection criteria. With a non-zero f°, and small CP violation in 
D° — D° mixing, the asymmetry will be diluted making the deviation from the SM even 
more severe. 



II. CONSTRAINTS ON NEW PHYSICS PARAMETERS 

The required value for is much larger than the SM prediction. Attributing the observed 
excess to a contribution from B s — B s mixing, one needs to explain what type of new physics 
can produce the value (a^avc ~ — (12. 7 ± 5.0) x 10~ 3 . Theoretically, in terms of the mixing 
parameters in the B s — B s system, to a very good approximation a s sX is given by 

a ^P#^W^ (497±0 ' 94)xir W- (6) 

where A s and the phase <p s are defined by, M s 12 = M S 12 - 5M + M S 12,JVP = M} 2 ' SM A S = 
\M} 2 ' SM \ | A s \e l ^ s . We adopt the phase convention that T^ 2 is real. 

bmce the SM prediction for AM S w 2\M] 2 \ agrees with data well, |A S | is only allowed 
to vary in a limited region, 0.92 ± 0.32 fixed by the experimental value [s| AM S = (17.77 ± 
0.12)ps~ 1 and the SM prediction Q (19.30±6.74)ps~ 1 . In order to reproduce the DO result, 
it would seem naively that sin <p s would have to exceed the physical range | sin <p s \ < 1, as one 
needs sin0 s = —2.56 ± 1.16. This situation calls for a more careful analysis considering <j) s 
and |A S | simultaneously. Here we wish to point out that Eq. ([6]) is only an approximation, 
and we now review the derivation of the exact formula for the asymmetry. 
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Denoting the element in the Hamiltonian responsible for B s — B s mixing by M s 12 — iT 12 /2, 
and working in a basis where — T] 2 = |rj. 2 | is positive real, one can write the same element 
as \Ml 2 \e i<t>a + i\T 12 \/2 (this is equivalent to defining (f> s = arg(-M] 2 /V 12 )). We have 



AM, + i/ST s /2 = 2^(M} 2 - iT 12 /2)(M} 2 * - iT 12 * /2) , 

AM S AT S = 4|M s 12 ||rf| cos0 s , AM 2 - Ar 2 /4 = 4(|M S 12 | 2 - |rf |/4) . 

Note that the above definitions of AM S and Ar s are the same as those in Ref. 
Further defining w s = 2\M} 2 \/T S , and z s = \Tl 2 \/T s , we have 



w 2 s - z 2 = x 2 - y 2 , w s z s cos 4> s = x s y s 



(7) 



(8) 



which lead to 



4:X 2 y 2 
cos 2 



ixi 



,2\2 



x« 



[Xt 



V. 



2\2 



Ax 2 y 2 



cos^ 
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The asymmetry in terms of w s , 2 S and sin S is given by [9 ] 

sin S 



2w s z s sin i 
w 2 + z 2 



(xj-yi) 2 cog 2 



sm i 



1 



(l-(Ar a /2AA/ s ) 2 



4(Ar s /2AM s ) 2 COS ^ s 

The fact that i« s >> z s allows one to write an approximate formula 



'si 



2z„ 



U2| 



sm i 



sm i 



w s \Ml 2 \ 

This is the formula that is often used and given in Eq. ([6]). 
A careful expansion of a s sl in terms of x s and y s , reveals that 



|^tan0 s 



sm i 



AT, 
AM S 

AT, 
AJ\i s 



<< COS ' 
>> COS ' 



(9) 



(10) 



(11) 



In the SM, the phase in M 



12.SM 



(12) 



is |2[ 0.0041 ± 0.0014 which is too small to play a 



Note that the asymmetry can be as large as order one for fixed AM, and Ar s . 



substantial role in explaining the large asymmetry observed. New physics beyond the SM 
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may induce large CP violating phases and also change the magnitudes for both M] 2 and 

r 12 . 

s 

We now consider the constraints on the new physics contribution to M] 2 ' NP = 
|^i2,ATP| e i0jvp assuming that there is no alteration to the SM prediction for T] 2 . Defin- 
ing R = \M} 2 ' NP \/\M} 2 ' SM \ and neglecting the small phase in M] 2 ' SM , we can solve for 



sine/), and A 



sine/), .Rsini 



>NP 



|A S | (I + 2R cos (j) NP + R 2 ) ' 

|A S | = (1 + 2R cos (j) NP + R 2 ) 1/2 (13) 

Whether there are physical solutions for R and 4>np, should be analyzed using the above 
equations with the constraints for sin0 s /|A s | and |A S | in the ranges (—2.56 ± 1.16) and 
(0.92 ±0.32), respectively. 

In Fig.l we show the ranges for |A S | and R for a given value of 8 = sin</> s /|A s | = 
a^/4.97 x 10~ 3 (or = 4.97 x 10 -3 <5) as a function of sin 4>np- Since DO data requires that 
the asymmetry a s sl be negative, this restricts sin cf>NP to be negative too. One can then see 
in which quadrant should <Pnp be in order to reproduce the data. To obtain a large size for 
5, a lower value of |A S | is preferred, therefore the solution with cos<p s < is preferred. 

From the figures, we see that it is not possible for 5 to get down to the central value 
—2.56 as this would require a value for |A S | below its one sigma lower bound. We have 
checked that in order to have solution within the one a region of |A S |, 5 can at most go 
down to -1.6. To have 5 reach the DO central value it is also necessary to modify r* 2 . The 
analysis, in general, will now be different 9J. However, since |M S 12 | is much larger than |r* 2 |, 
the change of Tj 2 needed can be easily accommodated in Eq. (jBJ) by multiplying by a factor 
9 = |rJ 2 |/|r}2 ,5M | and modifying the phase <p s to include the contribution from T^ 2 . The 
value for a s sl is then scaled by a factor 9. The central value of the DO asymmetry is then 
obtained with 9 around 1.6, which is allowed by the experimental data on Ar s . It should 
be noted that the usually quoted value for T 12 is from short distance SM contribution 
there may be long distance contributions which modify the value. A larger than SM short 
distance contribution to T^ 2 is still a possibility. However, it is difficult to reliably calculate 
the long distance contribution. It is also possible that a large T] 2 is due to new physics 
beyond the SM 
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When we go beyond the SM, new contributions are in general not known. It is therefore 
desirable to use experimentally measurable quantities as much as possible. The exact formula 
in Eq. (TTOT) allows one to predict using the measured values x s and y s , and a theoretically 
unknown phase <j) s . Taking the experimental values j^] x s = 26.2±0.5 and y s = 0.046±0.027, 
one can ask what theoretical values for |M S 12 |, |r* 2 | and sin0 s are needed. We show the 
results in FigfSJ 

The central value of dimuon asymmetry can be reproduced with sin <p s around - 0.96, and 
z s around 0.16. This implies that Tj 2 is a factor of 3 larger than the usual SM short distance 
contribution. 



III. IMPLICATIONS FOR BEYOND SM MODEL PARAMETERS 



In order to reproduce the anomalously large dimuon asymmetry observed by DO, new 
physics beyond the SM is needed. There are different ways in which the dimuon asymmetry 
in B s — B s can be affected. For example Z' models with tree level FCNC jol, 24] can 



easily induce a large modification for M 12 , but have little effect on V 12 . R-parity violating 
interactions in SUSY models can also have tree level FCNC. Besides modifying M 12 , they 



various experimental data 



23 



can induce modifications in T] 2 11 1. However, taking into account new constraints from 



3H26|. it is difficult to generate sizeable modifications in 
T 12 1 . In the following we discuss these two types of models and their contributions to M s 12 
and Tj 2 , taking into account constraints from other possible data, in some detail. 





-1.0 -0.8 -0.6 -0.4 -0.2 0.0 



FIG. 1: R and |A S | as functions of sin^jvp for representative values of 5. All curves use the central 
value for M 12 > SM . 



We than C. Bauer and N. Dunn for bringing these new constraints to our attention. 
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Z' model with FCNC 



In general a Z' can couple to quarks with FCNC of the form 24j 



C 



2cw 



(14) 



In Z' models, the new contributions to M\ 2 happen at tree-level, whereas the modification 
of r^ 2 only occurs at the loop level and is therefore a smaller effect. We will concentrate 
on the tree level Z' exchange contribution to M] 2 with constraints on the parameters from 
other processes 23l-l25l|. We begin from the known contribution to M12 from Z' exchange 



2i, 



M. 



z , Gp m 2 z 6/ 23 1 , 2 M rj / 2 , l2 

- ^rrfVz' ^JB s M B s ^B s {a sb + b sb 
+ %f /23 ^a sfe 6 s6 (2e - 3) + ^(%, 3/23 



Vz^ /23 )-J a sbb s b(l - 6e)). 



(15) 



where r\z> = a s (m^)/a s (/i) is a QCD correction resulting from the running scale from mz 1 
to \i — m^. Bb s is the ratio between the matrix element < Bslsj^^bs^nJ^Bs > and its 
value in factorization. Similarly, e is defined as e = (Blr/ BB a )(rn 2 B J {m s + rrib) 2 ) where Blr 
is the ratio between the matrix element < B s \s^^bs^j^b\B s > and its value in factorization. 
When needed, we will use e = 1 in our numerical results. 

Using the central value from lattice calculation (2?]] fB s \]P>b s — 270 MeV, we obtain 

2 



1 + 6.8 x 1(T 



mz 



x 



« 2 6 + b lb + Vz'^^sbbsb^e - 3) + jj(%? /23 



Vz 



30/23 



-a sb b sb (l - 6e) ) . (16) 




1.0 -0.9 -0.8 -0.7 -0.6 -0.5 
Sin I, 




FIG. 2: From left to right oA, w s = 2\M} \/T 8 and z s = |Tj |/r s as functions of sm.cj) s respectively. 
In all cases we use the experimental central values x s = 26.2 and y s = 0.046. 
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To make sure that there are solutions for the required value of determined from DO 
data, we take with b s b = for illustration. In this case we have 

R = 6.8 x 10 4 (—) \a sb \ 2 , <j) NP = 2arg{a sb ) . (17) 

As we have discussed earlier, within the one a allowed region for |A S | it is not possible 
to obtain the DO central value for a s A . To illustrate the range of a\ x that can be obtained 
with this model, we consider a few specific values of 5 = sin0 s /|A s | and restrict |A S | to be 
within its one a allowed region. 

A solution with an asymmetry within one a of the DO result requires 5 to be less than -1.4. 
For illustration we take 5 = —1.5 (corresponding to in the range of (—0.78 ~ —6.0) x 10~ 3 ). 
With 5 = —1.5, sin 0atp is restricted to be in the range —0.39 ~ —0.59 and the corresponding 
range for R is 1.37 ~ 0.91, as can be seen in Figured! If the large asymmetry observed by 
DO is confirmed, the new physics parameters R and sin <p^p need to be in the above ranges. 
However, if it turns out that the asymmetry is smaller, one needs to use a smaller 5, in which 
case solutions are much easier to obtain. We will concentrate on the case with 5 = —1.5. 
Applying the above ranges of parameters to Eq. (fT7"|) , we then have 



7T 7T 

arg{a s b) is in the range: (— + 0.19 ~ — + 0.31) , 

^\a sb \ is in the range: (0.0036 ~ 0.0044) . (18) 
m 7 



This range is comparable to other constrains on FCNC from non-universal Z' models 23], 
and in particular admits solutions with small |a s &|, say of O(10 -2 ). We note that Z' couplings 
to bs quarks has very weak constraints from the decays B s — > fi + fi~ and Bd — > Kl + l~ which 
are much lower than values used above. In this case the Z' mass can be in the several 
hundreds of GeV, a region that can be studied at the LHC. We note that models with 
natural suppression of flavor changing couplings exist in the literature 251 ] . 

The case with a sb = is the same as the case discussed above, but with a sb replaced by 
b sb . If both a sb and b sb are non-zero, the analysis is more complicated. For example, for the 
special case with a s b = b sb the contribution from Z' exchange to R is reduced by a factor 
of 0.6. This translates into the coupling, a s b = b s b, being enhanced by a factor 1.3 and the 
phase range being the same as in the b s b = case. 
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B. SUSY models with R-parity violation 

We now discuss an example which can modify T] 2 , an R-parity violating (RPV) interac 



tion in SUSY models. However, we find that existing constraints 
levels below what is required to explain the DO data. 



will limit the effect to 



There are three types of R- Parity violating terms [28[: X i jk/2L l L L :> L E R k , X' i - k L l L Q :1 L D'^, and 



\'^ k /2U R D c R D R . Here i, j and k are the generation indices: L L , Q L , E R , D R and U R are the 
left handed lepton, quark, right handed lepton, down-quark and up-quark fields, respective. 
ijj c indicates the charge conjugated field of ip The contributions of these interactions to ArJ 2 
have been studied in detail ll|. It is found that these interactions can induce a non-zero 
M 12 at tree level. There are also couplings that can induce a non-zero V 12 at one loop level 
without tree level contributions to M\ 2 . 

The contributions to z s can be grouped into several categories according to particles 
exchanged For A' couplings, these contributions are: z s (SM — RPV) from interference 
between SM and R-parity violating interactions; z s (RPV — RPV,v) from exchanges of 
neutrinos and down- type squarks; z s (RPV — RPV, I) from exchanges of charged leptons and 
up- type squarks; z s (RPV — RPV,u) from exchanges of light up- type quarks and charged 
sleptons; and z s (RPV — RPV, d) from exchanges of light down-type quarks and sneutrinos. 
For A" couplings, the contributions are: z s (SM — RPV) from interference between SM and 
R-parity violating interactions; z s (RPV — RPV,u) from exchanges of up-type light quarks 
and down-type squarks; and z s (RPV — RPV,d) from exchanges of down- type light quarks 
and up-type squarks. 

Although the couplings involved do not contribute to B s — B s mixing at tree level, they 
contribute to various B decays, such as b — > 57 and B — > MM (where M is a light pseu- 
doscalar meson). After applying these constraints, the most likely large contributions are 
given by 

(lOOGeV) 2 



z s (SM - RPV) = -13(A; 22 A:* 23 + A; 22 A 



/* 



z s (RPV - RPV, d) = -233 x 28 x K^mK^f ^f^ , (19) 



z s (SM - RPV) = -2.9A 221 A 2 '* : 



(lOOGe^) 



2 



221^231 m 2 
"l 
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The constraints on these couplings that arise from B decays, assuming sparticles with 
mass 100 GeV, are Q, |A- 32 Ag 2 | ~ |A^ 23 A£ 2 | < 2.3 x 10" 3 , |A^ 13 Af 22 | < 2.48 x 10" 3 , and 
I "^221^231 1 < 10~ 2 . These numbers constrain the three contributions in Eq. [19]to be less than 
0.06, 0.03 and 0.03, respectively. All of them are smaller by factors of 3-5 than the z s ~ 0.16 
required to explain the DO asymmetry. We cannot rule out the possibility that all these 
contributions (plus the SM) add up constructively to reach the value required by DO data, 
but at this stage it is fair to conclude that this is not a favored explanation. 



IV. TIME DEPENDENT CP VIOLATION WITH A NON-ZERO Ar 



Finally, we comment on the CP asymmetry A^cp which can be measured by studying 
the time dependent B — > l + uX and B — > l~vX decay rate difference, 

^ t AfCoa(AMsQ + S f sm(AM 8 t) 
AtCP = 26 i + e Ar s *_ A Ar (1 _ e Ar sf) > ( 20 ) 



where / is not a CP eigenstate. The notation follows Ref. 24] 



A _ W)! 2 - W)! 2 9 Im{{q B J PBa )A{f)A*{f)) 
f \A(fW + \A(f)r f ~ \A(fW + \A(fW ' 

Af = 2 g%#M^W :Q1 : \M 2 + \Sf\ 2 + \Af? = 1. (21) 



In this equation, A(f) and A(f) are time dependent decay amplitudes for B s and B s decay 
into states / and / in terms of the B s mixing parameters 





I M} 2 * 




Pb s 


V Ml 2 





(22) 



Assuming that CP violation in A and A is small, \A\ = \A\, one obtains 



at^ sin0 s sin(AM s t) 
ATCP ~ 2e l + e Ar.t- co8^(l -e^" (23) 



If Ar g = 0, which is a very good approximation for B^ decays, it is not possible to study 
the quantity A^ r in totality. The time dependence is a simple sine function of time and one 
cannot check the Ar effect. Therefore this CP violating observable is special to B s — B s 
system because Ar s is not zero. It offers a possibility to study the exponential factor in 
the time dependence. In the SM the phase <p s is small j^, 0.0041 ± 0.0014, resulting in a 
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very small Atcp- If the DO result is confirmed however, a large phase is allowed as we saw 
earlier and a large A TC p is possible. In FigQ, we show cltcp — A TC p(AT s ) — A TC p(0) as a 
function of t. We have chosen the value sin0 s = —0.5 with the phase in the third quadrant 
and the central value A s = 0.096 for illustration. We can see that at the few percent level, 
there are differences with respect to the A1T S = case, and such differences may be tested 
at the LHCb or at a B-factory, such as BELLE II with sufficient high energy to produce 
B S B S pairs. 

0.03 
0.02 
0.01 

b o.oo 

-0.01 
-0.02 
-0.03 
-0.04 

0.0 0.5 1.0 1.5 2.0 

t(ps) 

FIG. 3: cltcp vs. t(ps) with sm^ s = —0.5 and cos0 s < 0. 
V. SUMMARY 

The dimuon asymmetry reported by DO is much larger than the SM prediction although 
further experimental studies are needed to confirm this result. D° — D° mixing may con- 
taminate the final results and it is therefore important to carry out a detailed study with 
careful selection criteria for the dimuon events. If the DO result is confirmed, it represents 
clear evidence for new physics beyond the SM. 

If new physics only affects M s 12 , we can get an asymmetry within one a of the DO value, 
but it is not possible to reach the central value. Modification of T] 2 by new physics may 
then play an important role. We studied the consequences of the DO dimuon asymmetry on 
a Z' model with tree level FCNC and in a SUSY model with R-parity violating interactions. 
We find that the exchange of Z' can significantly modify M] 2 and bring the theoretical 
prediction to within one a of the DO allowed region. We showed that there are R-parity 
violating terms which can modify T] 2 but that these modifications are probably too small to 
account for the observed DO asymmetry. We also showed that the DO result implies a large 
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effect on time dependent CP violation in B s — B s mixing resulting in a detectable non-zero 
Ar s effect at the LHCb. 
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